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Asymmetric Heat Flow in Mesoscopic Magnetic System 
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The characteristics of heat flow in a coupled magnetic system are studied. The coupled 
system is composed of a gapped chain and a gapless chain. The system size is assumed to be 
quite small so that the mean free path is comparable to it. When the parameter set of the 
temperatures of reservoirs is exchanged, the characteristics of heat flow are studied with the 
Keldysh Green function technique. The asymmetry of current is found in the presence of a 
local equilibrium process at the contact between the magnetic systems. The present setup is 
realistic and such an effect will be observed in real experiments. We also discuss the simple 
phenomenological explanation to obtain the asymmetry. 
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1. Introduction 

Recently, quantum heat transport in magnetic materials has attracted much interest. In 
many magnetic materials, a magnetic contribution is dominant in heat transport in compari- 
son with phononic heat transport. Generally the thermal conductivity is sensitive to magnetic 
phase transitions. Therefore, measurements of heat flow provide information on magnetic 
properties.^ On the other hand, magnetic materials have also provided fundamental infor- 
mation on heat transport in nonequilibrium statistical mechanics. For example, which types 
of magnetic system show a diffusive transport in the thermodynamic limit? In an isotropic 
Heisenberg chain, the heat flow is a conserved quantity, and ballistic heat transport is theoret- 
ically expected.^ Generally integrable systems show a divergent conductivity.®'^ In Sr2Cu03 
and CuGeOs, which are well described by the isotropic Heisenberg chain, the ballistic feature 
of heat transport is actually confirmed by measuring an unusually large mean free path (~ 
500-1000 times the lattice const ant ).^"'^ Interestingly, the Heisenberg chain with extra pertur- 
bation terms also shows very large conductivity, which is shown for many materials such as 
the spin-Ladder systems (Sr,La,Ca)i4Cu2404i and La5Ca9Cu2404i.^"'' Although in real ma- 
terials, macroscopic heat transport in the thermodynamic limit cannot be ballistic because of 
the existence of spin-phonon interaction, impurities and so on, it is still of academic interest to 
study the thermodynamic behavior of magnetic heat transport. Therefore, both theoretically 
and experimentally, the properties of energy transport in the thermodynamic limit have been 
intensively studied. ^'^"^^ 
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In this paper, however, we do not study the thermodynamic hmit. We focus on mesoscopic- 
scale systems, where mean free path is comparable to system size. Therefore, phononic heat 
transport is negligible at very low temperatures, and magnetic heat transport is regarded 
as ballistic. We study a coupled magnetic systems and consider how the heat flow can be 
controlled. In this case, two magnetic materials are connected to each other, and it is assumed 
that one magnetic system is gapless, and the other is gapped. In this case, there exists no 
magnetic interaction between the two magnetic chains. Therefore, the heat is conveyed through 
phonons at the surface of contact between them. Here, we ask what should be expected for 
heat flow in this setup. 

A similar setup has been studied in classical anharmonic chains, where left and right chains 
have different phononic mode distributions. At the left and right edges, thermal reservoirs at 
different temperatures are attached. Terraneo et al. numerically showed that asymmetric heat 
flow is observed when the temperatures at both edges are exchanged.^''' This phenomenon 
is observed only for the regime far from equilibrium, and cannot be explained within the 
linear response theorem. Recently, it has been desired to demonstrate this phenomenon in real 
experiments.'^ In this paper, we propose a promising experimental setup to observe asymmetric 
heat flow using the above coupled magnetic system. It is emphasized that our setup is based 
on a mesoscopic quantum magnetic system so that we have a wide choice of magnetic material. 

Here, we theoretically analyze heat flow using the bosonization technique. The bosonized 
Hamiltonian is convenient because it guarantees ballistic energy flow, which is expected in a 
mesoscopic regime.^'' We adopt the Keldysh Green function technique for the boson operators 
to study the noncquilibriTini stationary state. We show that asymmetric heat flow is actually 
observed only when the thcrmalization of phonons occurs, which is normally expected in real 
experiments. 

This paper is organized as follows. In §2, we introduce the model, and analyze the station- 
ary state using the Keldysh Green function method in §3. In §4, another analysis to interpret 
the phenomenon is introduced. The summary and a brief discussion are given in §5. 

2. Coupled Magnetic Systems 

We consider the composite mesoscopic spin-^ system, which is composed of a gapless 
chain and a gapped chain. This attachment of two different mesoscopic magnetic systems 
means that heat transport between magnetic chains is possible by spin-phonon interaction at 
the surface of contact between them. Here, it is assumed that the mean free paths in both 
magnetic systems are comparable to system size. These setups are depicted in Fig. 1(a). We 
now use the following type of Hamiltonian to investigate the heat flow at the nonequilibrium 
stationary state: 

n = Hl + Hr + Hp + H' (1) 
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= l^'^fc^I^fc (2) 
k 

n' = ^ALfe4,i(5, + 6i)+5;AR,4,i(6, + 6t) (3) 

k k 
Here Ha {a = L and R) is the Hamiltonian of the left and right magnetic systems, and Hp 
expresses the phonon part. The operators sf^ ^ and s^^ ^ are the z-component of spins at the 
edges of the left and right magnetic systems, and they are assumed to interact with phonons 
localized at the contact between magnetic systems. These decompositions are schematically 
shown in Fig. 1(b). Here, the isotropic Heisenberg chain is assumed to be the gapless chain, 
but no concrete gapped chain is assumed at present. 

We treat mesoscopic-scale systems so that mean free path is comparable to material size 
in both magnetic materials, where heat is conveyed without scattering within the material. 
In this case, the effective Hamiltonian is well described by the bosonization form of the spin 
systems. Generally, bosonized Hamiltonians satisfy a ballistic feature of transport because 
the bosonized heat flow operator is commutable with the Hamiltonian.^^ Thus, the bosonized 
form of the spin system automatically satisfies the characteristics of the mesoscopic nature of 
the magnetic systems. We adopt the following bosonized Hamiltonian: 



u 
2 



poo |- -, 

/ dx n2 (x) + (a^$L(x))2 (4) 
Jxo 

Hn = - dx Ulix) + {dAix)f + go^iix) (5) 




^' = 2. 7^^-^«M + )^^(6fc + 6l) ) (6) 

a=L,R V"^ k 

The boson operators ^a{x) and Ilce'{x') satisfy the commutation relation [^oi{x),^a'{x')] = 
iSa,ce'S{x — x'). The operator -^dx^a{xo) is the z-component of spin at the position xq in the 
ath spin chain. In the interaction Hamiltonian H', the quadratic terms of bosons are added 
only to avoid an unbounded state. For simplicity of calculation, we adopt the bilinear form for 
the gapped Hamiltonian, which can be regarded as the approximated form of the sine-Gordon 
Hamiltonian Hr.^'^^ The sine-Gordon Hamiltonian can be mapped from many gapped mag- 
netic chains. The present gapped Hamiltonian gives the dispersion relation with the energy 
gap A = u^/go- We believe that this simplification does not cause different qualitative results 
because correct excitation spectrum and excitation eigenfunctions are not necessary for the 
phenomenon described in this paper. 

3. Heat Flow 

In this section, we calculate the heat flow at the steady state. As explained in the previ- 
ous section, the mesoscopic nature of current is represented by the ballistic heat flow in the 
bosonized Hamiltonian. The next task is the realization of the nonequilibrium steady state of 
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the system. For simplicity of calculation, we assume an infinite size for the left and right mag- 
netic systems, instead of connecting thermal reservoirs at the edges of the magnetic systems. 
By taking the infinite size, we can realize a nonequilibrium steady state as follows (see Fig. 
1(b)). We assume that both magnetic systems are initially separated and are in equilibrium 
with the temperatures Tl and Tr. Next, these systems are adiabatically connected with the 
phonon system, and a steady state is realized in the long time limit. The steady state can be 
analyzed by the Keldysh Green function technique. To this end, the following operators are 
defined 

= XI (^"'^ + ^ak) ' Ya = i^Xak (bak - Oi = R OT L (7) 

k k 

We further assume the simple form of coupling for spin-phonon interaction as Al^ = Ar^ = A^. 
We do not discriminate between the indices a of operators and Y^, and simply write X 
and Y, respectively. The heat flow is expressed by the Keldysh lesser Green function 



(8) 



Throughout this paper, we set H to be unity. Here Gyxi^) ^^'^ ^yx(^) Fourier 
transforms of the retarded and lesser Green functions: 

G'v^xM = -ieit-t'){[Yit),x{t')]) 

G<^{t,t') = -i{X{t')Y{t)) 

The functions ^{Cl) and g'^°g(f2) are the Fourier transforms of the nonperturbed Green 
function: 

glsit^t') = ~i{Bit')A{t))o 

9XbM = r^Q{±{t-t')){[A{t),B{t')])o 

Here (...)o denotes an average over a nonperturbed Hamiltonian. The nonperturbed Green 
functions are readily calculated as 



9xIyS^) 
9yIyS^) 
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T4,sgn(n)VJ22_A2 for |Q| > A, 

V^^^ for \n\<A, 



(9) 



(10) 
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where a is L or R, and Al = 0, and Ar = u^/qq. Here, we neglected the effects of the 
quadratic terms of the bosons in the Hamiltonian (6), which is not expected to have a large 
contribution to the stationary state. 

In realistic systems, it is natural to expect that the phonon system at the contact be- 
tween magnetic systems becomes a local equilibrium state at the steady state. To realize the 
thermalization of the phonon part, we adopt the so-called self-consistent reservoir (SCR) by 
connecting an extra infinite system, which can be regarded as the reservoir with some temper- 
ature Ts at the initial time. The temperature Ts is determined under the condition that heat 
flow vanishes between the reservoir and the phonon system at the stationary state as depicted 
in Fig. 2(a). Such an extra system induces the thermalization of phonons, and the degree of a 
local equilibrium state is controlled by the coupling strength between the SCR and the phonon 
system. Since this extra system is simply adopted only to induce the thermalization of the 
phonon system, we have many choices for the system. For simplicity, we now use the spin 
system for the SCR, which is the same form of because all Green functions are already 
known. The interaction Hamiltonian where the phonon couples with the SCR is 

WSCR-P = A ^-^a^$sCR(xo)+ + ' (12) 

where -^dx^sCRixo) is the z-component of spin at xq in the SCR, and A is the cou- 
pling strength between SCR and the phonon system. The strength of A controls the de- 
gree of local equilibrium. The spin system for the SCR is the same form as that T^l! 
WsCR = I /o°° dx [n|cR(x) + {dx^sCK{x)f \ . The Green functions and G^_x(^) 

in the formula (8) are calculated by the Keldysh relation^^ with the nonperturbed Green 
functions of the phonon system, two magnetic systems and the system of SCR. 

Let us now sec the properties of heat flow in these setups. Two cases are investigated, i.e., 
(IljIr) = (AT,0) and (Tl,Tr) = (0, AT). The heat flows in the former and the latter cases 
are written as Jl and Jr, respectively. In these two cases, the behavior of heat flow is inves- 
tigated for various coupling strengths A. Numerical calculation is sclf-consistently performed 
to satisfy the condition of SCR. Here, we took 1000 boson numbers for the phonon system 
with the ohmic spectral density. The energy gap is A = 3.0 in the unit of n. The temperature 
dependences of heat flow are shown in Fig. 2(b) for various coupling strength A. As shown in 
Fig. 2(b), in the case of A = 0, no asymmetry of heat flow is observed. However as A increases, 
asymmetry appears. For a larger A, the amplitude of heat flow decreases, which is due to the 
scattering effect of energy by the SCR. In the case of finite A, |Jr| is always larger than | Jl|. 

4. Extreme Model to understand Mechanisms 

As shown in the previous section, asymmetric current can occur in the presence of the 
thermalization of phonons. This means that the phonon system has some temperature de- 
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fined in the local equilibrium state. In this section, we consider a simple phenomenological 
explanation of the asymmetric heat flow regarding the phonon system as the phonon reservoir. 

To this end, the following two problems are considered, i) First, what amount of the 
stationary heat flow is obtained in one boson system when it is coupled to two different 
reservoirs with different temperatures ? Here the boson system is the main system of interest, 
ii) Second, the following situation is studied using the general formula of heat flow derived in 
the first problem. Suppose that the different two boson systems sandwich a phonon reservoir 
with some temperature, Ts. At the edges of the boson system, other phonon reservoirs are 
connected with the temperatures Tl and Tr (see the schematic explanation of the situation 
in the inset in the Fig. 3). What is the characteristic of heat flow, when the temperature Ts 
is determined under the condition that both stationary heat flows within the boson systems 
are equal ? This extreme situation actually explains the asymmetry of heat flow. 

We study the first problem by considering the dynamics for the Hamiltonian; 



n 



= ^+ H (^a-Il7a,fe(oa,fe + aL,fe) +J2^''"i,k^^'k, (13) 

a=l,2 \ k / a,k 

n = J^^kblbk, (14) 

k 
k 

where Ti is the boson system of interest and the operators aa,k and o)^ ^ are the annihilation 
and creation phonon operators in the phonon reservoirs. The operators Xi and X2 are the 
system's operators attached to the phonon reservoirs. We assume that the phonon reservoirs 
have spectral densities /a(^) and temperatures Tq,. 

By standard projection operator technique,^'' we obtain the Redfield-type master equation 
for a reduced density operator^^'^^ up to the second order of v 

^ = -i[n,p{t)]-v^C,p{t)-v^C2p{t) (16) 

C^p{t) = [X„,i?,p(t)] + [X„i?„/>(t)]t (17) 
Ra = ^Mfe/a(^^it)na(i^fc)(e^«"'=6fe + 6i) (18) 

k 



Here, nc{^k) is the Bose distribution, nci{^k) = l/(e^"^'= — 1). When wc expand the stationary 



density matrix as pst = P^^ + '^^Pst^ + ' ' ' ; we obtain the zeroth order of density matrix as 



(0) „ e-'^^^'^uo^ok 1 
^st = , Bk = —\n 



(19) 



The heat flow operator from the first reservoir to system is derived from the continuity 
equation of energy 

j = u^nci. (20) 
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Thus, using the zeroth order stationary density matrix p^^\ we obtain the leading order of 
average current (j) = v^Tt ^HCip^^^^ . The average current (j) is readily calculated, 

(j) = y dsV{e)e 

where 'D{e) is the density of states at the energy e, and Pa{£) = /^^ — f^fe)- The 
function of Ha^k depends on detailed systems. Below we consider the special case with p\i.= 
1^2 for simplicity. This case is actually realized in some boson systems.^'' In this case, eq.(21) 
becomes simpler because of = /^2(£)- 

We are now in the position to discuss the second problem and clarify the essential mech- 
anism of asymmetric heat flow. Suppose that the different boson systems are connected via 
a phonon reservoir with temperature Ts, and are also attached to phonon reservoirs with 
different temperatures, Tl and Tr (see the inset in Fig. 3). We call the heat flow within the 
left and right boson systems as j(Tl,Ts) and j{Ts,TR), respectively. At the stationary state, 
it is demanded that both currents are equal, i.e., 

jin,Ts) = j{Ts,Tn), (22) 

from which the temperature Ts is self-consistently determined. Note that the heat flow is 
measured under these conditions. 

We now demonstrate that the asymmetry of current actually occurs when the energy 
densities of the left and right systems, X'l(£) and ^^{e) are different. We consider the simplest 
case of I'L(e) = doe and ^^{e) = doVe^ + A^. In this case, we calculate the energy flow Jl for 
the temperature set (TL,rR) = (Ar,0), and Jr for (TljTr) = (0, AT). In Fig. 3, we present 
the results of Jr and Jl under these conditions. Here, energy gap A = 3.0, and we used the 
same phonon reservoir with ohmic spectral density. As shown, we find that the asymmetry 
appears in the heat flow, and the qualitative results in the previous section are reproduced. 

5. Summary 

We discussed the possibility of observation of asymmetric heat flow in a coupled magnetic 
material. We particularly focused on mesoscopic-scale magnetic systems, where mean free path 
is comparable to system size. We found that an asymmetric current is observed as a result 
of the thermalization of phonons at the contact between magnetic materials. The present 
experimental setup is realistic, and we have many choices of realistic magnetic materials. The 
asymmetry of heat flow is a key ingredient for heat pumping. It will be important to discuss 
a realistic setup of heat pumping and the efHciency,^^ although we do not discuss it here. 

In a gapped spin chain, a phononic contribution may appear in realistic experiments. Ac- 
tually in CuGe03, which shows spin-Peierls transition, there exists a crossover from magnetic 
heat transport to phononic heat transport when the temperature decreases across the spin- 
Peierls temperature.^ The quantitative properties depend on spin-phonon interaction in the 



^l(£)Il(£)/7.2(£)/2(£) 

;Ui(e)/i(e) + ^2(e)/2(e) 



(n2(£) - ni(e)) , 



(21) 
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bulk materials. Thus, to realize spin heat transport in all temperature regimes, we must find 
a suitable material with very weak spin-phonon interaction. 

In the present study, we consider small magnetic systems, where the analysis becomes easy. 
In the case that mean free path is much smaller than material size, treatment is much more 
difficult than the present analysis because most gapped magnetic systems are nonintcgrable 
systems, and use of the bosonized Hamiltonian may not be justified. It will also be important 
to study such a region, because it is not clear whether observation of asymmetric heat flow 
requires the ballistic properties of current. To answer this question, real experiments would 
be very useful. 



8/13 



J. Phys. Soc. Jpn. Full Paper 

References 

1) A. V. Sologubcnko, K. Gioanno, H. R. Ott, A. Vietkine and A. Revcolevschi: Phys. Rev. B 62 
(2000) R6108. 

A. V. Sologubenko, K. Gioanno, H. R. Ott, A. Vietkine and A. Revcolevschi: Phys. Rev. B 64 
(2000) 054412. 

Y. Ando, J. Takeya, D. L. Sisson, S. G. Doettinger, I. Tanaka, R. S. Feigelson and Kapitulnik: 
Phys. Rev. B 58 (1998) R2913. 

M. Hofmann et al: Phys. Rev. Lett. 87 (2001) 047202-1. 

A.V. Sologubenko, K. Gioanno. H. R. Ott, U. Ammerahl, A. Revcolevschi, D. F. Brewer and A. 
L. Thomson: Physica B 284-288 (2000) 1595. 

A. V. Sologubenko, K. Gioanno, H. R. Ott, U. Ammerahl and A. Revcolevschi: Phys. Rev. Lett. 
84 (2000) 2714. 

K. Kudo, S. Ishikawa, T. Noji, T. Adachi, Y. Koike, K. Maki, S. Tsuji and K. Kumagai: J. Phys. 
Soc. Jpn. 70 (2001) 437. 

X. Zotos, F. Naef and P. Prelovsek: Phys. Rev. B 55 (1997) 11029. 
K. Saito, S. Takesue and S. Miyashita: Phys. Rev. E 54 (1996) 2404. 
J. V. Alvarez and C. Gros: Phys. Rev. Lett. 89 (2002) 156603. 

J. V. Alvarez and C. Gros: Phys. Rev. Lett. 88 (2002) 077203. 
S. Fujimoto and N. Kawakami: Phys. Rov. Lett. 90 (2003) 197202. 
E. Orignac, R. Chitra and R. Citro: Phys. Rev. B 67 (2003) 134426. 

E. Shimshoni, N. Andrei and A. Rosch: Phys. Rev. B 68 (2003) 104401. 
X. Zotos: Phys. Rev. Lett. 92 (2004) 067202. 

C. Hess P. Ribeiro, B. Buchner, H. ElHaes, G. Roth, U. Ammerahl and A. Revcolevschi: cond- 
mat/0506595. 

P. Jung, R. W. Helmes and A. Rosch: cond-mat/0509615. 

F. Heidrich-Meisner, A. Honcckcr, D. C. Cabra and W. Brenig: Phys. Rev. B 66, (2002) 140406(R). 
K. Saito: Phys. Rev. B 67 (2003) 064410. 

M. Terranco, M. Peyrard and G. Casati: Phys. Rev. Lett. 88 (2002) 094302. 

B. Li, L. Wang and G. Casati: Phys. Rev. Lett. 93 (2004) 184301. 

D. Segal and A. Nitzan: Phys. Rev. Lett. 94 (2005) 034301. 

L. V. Keldysh, Zh. Eksp. Theor. Fiz. 47 (1964) 1515 [Sov. Phys. JETP 20 (1965) 1018.]. 
A.O. Gogolin, A. A. Nersesyan and A.M. Tsvelik: Bosonization and Strongly Correlated Systems 
(Cambridge University Press, 1998). 

T. Nakano and H. Fukuyama: J. Phys. Soc. Jpn, 49 (1980) 1679. 
M. Bolsterh, M. Rich and W. M. Visscher: Phys. Rev. A 1 (1970) 1086. 

R. Kubo, M. Toda and N. Hashitsumc: Statistical Physics II (Springer- Verlag, New York, 1985). 
A. G. Rcdficld: Adv. Magn. Rcson. 1 (1965) 1. 

K. Saito, S. Takesue and S. Miyashita: Phys. Rev. E 61 (2000) 2397. 



For example, in case of the quantum harmonic chain whose Hamiltonian: H = J2n=i 



2m 



Y,n=o ^^^i^n+i - Xn)^ {xq = xn+1 = 0), Xi and X2 are xi and xn- The functions /ii,fe and iJ.2,k 
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1, 2, • • • , A^). These functions satisfy /zj = fc- 
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Fig. 1. (a): Experimental setup of coupled magnetic system. It is assumed that no magnetic interac- 
tions exist between the two magnetic systems. The heat can be conveyed by plionons at the surface 
of the contact, (b): Theoretical setup. At the contact between the magnetic systems, spin-phonon 
interaction exists. 

Fig.l, Keiji Saito 
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Fig. 2. (a) Schematic explanation of structure of self-consistent reservoir. Between the SCR and 
the phonon, there is no heat flow at the stationary state. Thermalization occurs at the phonon 
part with a flnite coupling strength, A. (b) | Jl| and | Jr| as functions of AT for various coupling 
strengths A of SCR. Here, the energy gap of the right spin chain is 3.0 in the unit of u. The unit 
of the heat flow is u'^/h. 



Fig.2, Keiji Saito 
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Fig. 3. Simple explanation of asymmetric heat flow using phonon reservoirs with temperatures Tl, 
Tg, and Tr. The temperature Tg is determined by the condition J(Tl,Ts) = J(Ts,Tr) 



Fig.3, Keiji Saito 
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